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Qv ■ In this paper we study the modification of the properties of the nucleon in 

CN , the nucleus within the quark-sohton model. This is a covariant, dynamical 

model, which provides a non-linear representation of the spontaneously bro- 
^ \ ken SU{2)l x SU{2)r symmetry of QCD. The effects of the nuclear medium 

1^ ' are accounted for by using the Wigner-Seitz approximation and therefore re- 

(^ ■ ducing the complex many-body problem to a simpler single-particle problem. 

^— ^ , We find a minimum in the binding energy at finite density, a change in the 

Q"^ , isoscalar nucleon radius and a reduction of the in-medium pion decay con- 

^* \ stant. The latter is consistent with a partial restoration of chiral symmetry 

at finite density, which is predicted by other models. 
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I. INTRODUCTION 

In this paper we want to address the possibihty that the nucleon properties be modified 
in the nuclear medium. In a conventional nuclear physics approach, where nucleons and 
mesons are elementary degrees of freedom, such question cannot be answered consistently 
and any modification of the nucleon properties has to be put in by hand. In order to 
address the problem, one has to consider models in which the substructure of the nucleon 
is not neglected (see, e.g, Refs. f!]^), and properly implement these models in order to 
account for the presence of a medium. 

In the present work we have considered a chiral model of the nucleon, which has been de- 
veloped by Diakonov et al. [^,^ on the basis of the instanton picture of the QCD vacuum. It 
provides a low-energy approximation to QCD that incorporates a non-linear representation 
of the spontaneously broken chiral symmetry. In this framework pions emerge as Goldstone 
bosons, dynamically generated by the Dirac sea. Vacuum fiuctuations (quark loops) are 
described by an effective action that yields the pion kinetic term, — which is already in- 
cluded at the classical level in the Lagrangian of other chiral models, — and higher order 
non-local contributions. The model is also Lorentz covariant and has essentially only one 
free parameter (apart from the regularization scale), namely the constituent quark mass. Al- 
though the latter should in principle be momentum dependent, in practice a constant value 
is usually chosen, which better reproduces the phenomenological properties. The model has 
been successfully applied to the description of a variety of nucleon properties [Q-|^ . 

Solving the many-body problem is already a formidable task in a conventional nuclear 
physics approach, the more so, of course, when one deals with extended objects. Early 
attempts treated nuclear matter as a crystal, letting the particles sit on a regular lattice 
UsHTTI . However, one has still to face serious computational difficulties in properly imposing 
the Bloch boundary conditions and moreover, nuclear matter does not show long-range 
crystalline order. These facts prompt the application of an approximation first introduced 



by Wigner and Seitz |T^, in which the effect of the surrounding matter on each particle 
is accounted in an average manner, by enclosing it in a spherically symmetric cell: This 
technique does not depend on any particular structure of the lattice and it is particularly 
suitable for nuclear matter, which may be pictured more as a fluid than a crystal. In this way, 
a complex many-body problem is reduced to a single particle problem, where the effects of 
the nuclear medium enter only through average boundary conditions. The long-range order 
implied by imposing periodic boundary conditions gives rise to a band structure of the energy 
levels and one has to choose suitable boundary conditions for the lowest and highest energy 
levels. 

The Wigner-Seitz approximation to the treatment of soliton matter has already been ap- 
plied using other models of the nucleon structure: The Skyrme model 0,0], non-topological 



soliton models |T^-|T^, the hybrid soliton model |TB|- p(]|JT7|] and the global color model 
2T|^3[. In these models the boundary conditions for the spherically symmetric bottom 



level of the band[| are an extension to a Dirac spinor of the non-relativistic requirement of 



^ For the top level the problem is complicated by the lack of spherical symmetry (see Ref. [|17|). 



having a flat wave function. Moreover, in tlie cliiral soliton models |T7|-pO| the requirement 
of unit topological number inside the cell is also taken into account through the boundary 
condition on the chiral angle. 

In this work we explore the sensitivity of the calculation to the choice of different bound- 
ary conditions, by also imposing the requirement of flatness on the chiral angle, as in 



Refs. [|,|ll]. We also show that it does not make sense to discuss the band structure of 
the nuclear system without accounting for the spurious contribution to the energy stem- 
ming from the center-of-mass motion of the bags, since the corrections turn out to be larger 
than the band width and strongly dependent on the boundary condition for the level. 

Another important difference of the present work from the above mentioned calculations 
is the fact that here we also account for contributions coming from the Dirac sea. Since the 
chiral field is a mean field, it means that we include one-loop quark fluctuations. The latter 
generate the kinetic pion term, — which is included by hand in the model Lagrangians of 
Refs. p!7|-pn|, — and a (attractive) non-local contribution, which is known to be important 
for the calculation of free nucleon properties @-0. Dirac fluctuations are usually neglected 
in the in-medium calculations, although there is no reason to expect they are negligible or 
at least independent of the density. 

The paper is organized as follows: In Sections p.lA| and |111J| we briefly discuss the main 



features of the quark-soliton model and the approximations employed in its implementation; 
in Sections |II CHII Hi we introduce the Wigner-Seitz approximation, the appropriate bound- 



ary conditions for the fields and show how a few observables can be calculated in this model. 
A new orthonormal and complete basis in the elementary cell is also obtained, in which 
physical quantities, such as the vacuum energy, are expressed. In Section |T1| we present the 



numerical results, obtained by solving the equations of motion. Finally, in Section |^ we 
draw our conclusions and discuss possible future developments of the model. 

II. NUCLEONIC AND NUCLEAR MODELS 

A. Chiral quark-soliton model 



The chiral soliton model 0-y], which provides a non-linear representation of the SU{2)l x 
SU{2)r symmetry of QCD, is based on the lagrangian 

C=^[ip~MU5{x)]ip, (2.1) 

where ip represents the quark fields, carrying color, flavor and Dirac indices, while f/5 is a 
chiral field defined as 

U,{x) = i±^ U\x) + i^ U{x) (2.2) 

U{x) = exp [iT ■ 0{x)] . (2.3) 

The large (= 350 MeV) dynamical quark mass M which appears in Eq. (p.l|) is the result 
of the spontaneous breakdown of chiral symmetry, which also accounts for the appearance 
of massless Nambu-Goldstone pions. In this model the nucleon emerges as a bound state 
of Nc quarks in a color singlet state, kept together by the chiral mean field. Note that no 



explicit kinetic energy term for the pion is present in ( |2.1|) : Actually, the ip and U fields 
are not independent and the latter is in the end interpreted as a composite field in the 
quark- ant iquark channel. 

Introducing the familiar hedgehog shape of the soliton, 

f/(x) =exp{iT-f 0(r)}, (2.4) 

the quark hamiltonian reads 

H = -ia ■ V + PM [cos 9{r) - i-fr^r ■ r sin 9{r)] . (2.5) 



However, solutions to the classical field equations derived from ( |2.5| ) do not describe nucleon 
and A states, since angular momentum and isospin do not commute with H: One defines 
the so called grand spin, K = J + t/2, for which [H, K] = 0, and the quark wave function 
can be written as 



where ^ is the grand spin state fulfilhng 

(cr + r) e = (2.7) 

and the normalization is 

Att / dr r^ ^^7^ = 1. (2.8) 

Jo 

Good spin and isospin quantum numbers may be obtained in the end in a semiclassical 
approximation by quantizing the adiabatic rotational motion in isospin space |24,2|. How- 



ever, in the present paper, for the sake of simplicity, we limit ourselves to consider soliton 
matter, leaving the projection of spin-isospin quantum numbers for future work. 

The classical solutions are found self-consistently by solving the equations obtained by 
minimization of the total energy 

Etot [V, ?, 0] = N,E,,i [ij, ^, 9] + Evac [9] , (2.9) 

where E'vai = (V^I-f^lV') and E'vac are the valence and vacuum part of the energy, respectively. 
The vacuum energy E^^c incorporates, at the mean field level, one-quark-loop contributions 
and, formally, can be evaluated through the effective action, which is obtained by considering 
the following path integral over the quark fields 

exp {iSesp]} = f[d^][dilj] exp (i j d^x ^{x) [ip - MU^i^x)] %l){x) 

= [det {ip - MU5{x))f' . (2.10) 

The latter can be easily cast in a more suitable form by means of simple algebraic manipu- 
lations 



5eff[f/] = -^iV,trlog 






(2.11) 



where the trace is over Dirac and flavor indices. Despite its apparent simphcity, ( |2.11| ) is 
actually a complicate nonlocal object. 

Although a local derivative expansion of course is possible, it is of little practical use 
in this case, since the soliton fleld turns out to vary signiflcantly over the relevant distance 
scale M~^ and no stable solutions are found for expansions up to sixth order in derivatives 
p5| . Kahana and Ripka [^] on the other hand have developed a numerical algorithm to 



directly evaluate vacuum polarization contributions to soliton observables. This technique 
has been extended in Ref. @] to the calculation of nucleon observables, that is after collective 
quantization has been applied to project out states of deflnite spin and isospin. 

B. Effective action up to second order in the soliton field 



Another path, to which we shall adhere in the following, has been followed in Refs. PJ27[ 
by expanding (|2.11|) up to second order in M^pU^), obtaining 

Si^[U] ^ '-N^tr (x|^^ m^U,] ^^ m^U,] \x) 



-jN.ti f d^xd^x'V{x)K{x,x')V{x'), (2.12) 



where 



V{x) = iMpU^ 
K{x, x) = -tG{x, x)G{x', x) (2. 13) 

{n + M^)G{x,x') =6\x-x'). 

Note that in the standard derivative expansion the second order action would read 



In contrast to (|2.14| ) , Eq. ( p.l2|) does not assume a slowly varying soliton fleld and gives 



rise to non-local contributions. Furthermore, one can see that it gives a good approximation 
both for small and large momenta, thus providing an interpolation formula between these 
regimes [|,|1. 

Specializing to static fleld conflgurations, introducing E^^ = —S^n/Jdx'^ and going to 
momentum space, one gets 



where 



O J\/f2 

tT[V{q)V{-q)] = -—q^[Mq)<Po{q) + Mq)M-Q)], (2-16) 

•In 



with 



00 (g) = 47r/^ / drr'^jo{qr)[cos9{r) - 1] 



and 



[q) = iQiA-Rf^ / drr ji{qr) sin 6'(r) = iq^q), 
Jo 



K{q) = dre'''''K{r) 



(2.17a) 
(2.17b) 



(2.18) 



K{r) 



I /.oo .oo 7,- C\(k + k') /2] .,,,.,,, , 
dkk^ I dk'k'^^ ^^ ^r-^-^Joikr)jo{k'r). 



^vr^ Jo 







EhEu Eu + Ey 



In ( |2.18| ) Ek = y/k'^ + M"^, whereas C{k) is a regulating function, which will be discussed 
later. 

Although here and in the following we display, for convenience, formulae using a mo- 
mentum cut-off regularization scheme, we shall employ also the Pauli-Villars regularization, 
also to be discussed later. 

EI^^ must contain the meson kinetic energy contribution, which can be seen ||2^ to cor- 
respond to keeping only the g = term in an expansion of K{q) in (|2.15|) . This requirement 
fixes the normalization of K in such a way that 



1 r°° 

m^-,l dkk'- 



'C{k) fl 



El iN^M^ ■ 



(2.19) 



Then one finds E^^ = E^'"" + E^'^\ where 

\ /"OO 



E 



kin 



2Tifl drW'^e'^ + 2sm^ei 
Jo 



with 6' = dO/dr, and 



dqq"^ 



K{q) 
K{0) 



[<Pliq) + ^'m- 



47r2 Jo 

The propagator K{q) can be brought (see Appendix ^ into the form 

1 



Kiq) 



in'^q Jo 



■*Ar,,'*'c^i(*+*')/2i 



Ek J\k-q\ Ey Ek + Ek 



(2.20) 



(2.21) 



(2.22) 



C. Nuclear matter in the Wigner-Seitz approximation 

In order to describe nuclear matter we shall employ, as anticipated in the Introduction, 



the Wigner-Seitz (WS) approximation |[T^, which amounts to enclose the fields in a spheri- 
cally symmetric cell of radius R, imposing suitable boundary conditions. Before discussing 
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our choice of boundary conditions, let us describe the evaluation of the vacuum energy in 
the WS cell. 

We have first to find an orthonormal and complete basis of functions inside the elementary 
cell. We have chosen a spherical basis, in which the radial dependence is expressed through 
spherical Bessel functions, which have vanishing derivative at the boundary. This is the 
most useful basis in which to perform the calculation, when flat (zero derivative) boundary 
conditions for the fields are invoked. More details on the basis are given in the Appendix 
0. All the quantities involved in the calculation of the vacuum energy turn out to converge 
quickly using this basis. They also converge rapidly when zero boundary conditions are 
employed for the fields. 

Starting again from the static limit of (p.l2|) and (|2.13| ), one can introduce the Bessel 



transform of K{r, r') as 

K{r,r') =Y.Y: YUr)Y;m{r)Pa,ir)Pa'Xr')Kiiai,a[), (2.23) 



Im 



aia 



where 



Pa,{r) = f^aMc^ir/R), %t=^ = 0, (2.24) 



dr 
Kn, being a normalization constant. Inserting (|2.23|) into ( p.l2|) , one finds 



"a; 



-^ws — -^ws + -^WS' (2.25) 

where, as in the previous subsection, the kinetic energy contribution, — the one stemming 
from the local part of K{r,r') = 6{r — r')KQ{r) + ..., — has been separated. Indeed, one 
has 

rR 



with 



E^l = SttN.M^ f drKo{r)\r^e'\r) + 2sm^e{r)], (2.26) 

Jo 



^o(^)-i:i:^%^pU^) (2.27) 



and E^jR = ^{ai/Rf + M\ 

In the limit R —>■ oo one can check that 



Koir) ^^ Tdk e^ ^ -S—. (2.28) 



57r2 Jo El ANcM"^ ' 

As in the free space case, this fixes the normalization of -ft'o(^) in such a way that 

Then, one can write 



E. 



kin 



« , K.ir 



^« - '"^'^ L "^''iwy 



„2/]/2 



r'e'"{r) + 2sm'e{r)]. 



As we shall discuss below, one can view 



Jtt,WS — Jn 



2^o(r) 



Ko{0) 



(2.30) 



(2.31) 



as an in-medium, r-dependent pion decay "constant". 

On the other hand, the non-local contribution i^^S; — the in-medium extension of ( p.21D 
— can be cast into the following form: 



E. 



(2) 



ws 



— 7§ — Yl {i^foiao)AKo{ao, ao)/o(ao) + fi{ai)AKi{ai, a[)fi{a[] 

Jit aa' 

2 

+ o/2("2)Air2(a2, a2)/2(a2)}, 



where 



/o(«o) = Mi/2 / rfrrVoW 



R 



cos9{r)9\r) + 2 



sin 6'(r) 



R 



/i(ai) = M^/^ / drr^p^^{r) sme{r)9'{r) 



Jo 



R 



cose{r)e\r) 



smy r 



and 



°^ LL' \ 



L L' I 



KiLL'[(yi,(yi) 



E 



n 



E^ 



rR 



0^ 



f^lLL'{.(^ha'l) 



(2.32) 

(2.33a) 
(2.33b) 
(2.33c) 

(2.34a) 



i-ti 

^ u 



2^2 



^(a/,aL,aL/)^(a;,aL,aL/) (2.34b) 

(2.34c) 



One should now notice that the straightforward application of ( ^.11| ) to the Wigner-Seitz 
cell is incomplete, because it does not account for the Casimir energy intrinsically connected 
with the change of topology, which is present even in the absence of background fields. As 
a matter of fact one should write: 



^effft/] = --iV.tr {log [n + M^ + iMpU^{x] 

'n + M'^ + iM^U^ix) 



ws 



l°g[° + ^1.ee} 



--A^ctr log 



n + M2 



I A QCasimir 



Jws 



where 



A oCasimir 



^eff 



^iV.tr {log [n + M']^^ - log [n + M'lJ 



(2.35) 



(2.36) 



Note that the first term in ( p.35| ) is just the apphcation of Eq. (|2.11| ) to the Wigner-Seitz 
case; on the other hand, the second term is the genuine Casimir energy due to the change 
of configuration space. By performing the intermediate algebra one obtains 




^^casimir ^ 4^^y ^ EE(2^ + 1)\/^' + ^ - ii^EE(2^ + 1)V^' + ^ • (2-37) 

The Casimir energy, which is obtained by dividing this expression by the time — T, has 
now the expected form: It is the difference between the zero-point energy in the finite 
volume, obtained by filling all the negative energy orbitals in the Dirac sea, and the same 
expression in free space. As it stands, however, Eq. ( p.3?l ) is badly divergent and needs to be 



regularized. Such a task has been recently carried out for massive fermions, in the context 
of the MIT bag model, in Ref. [^] by means of the zeta function regularization technique. 
In the present work we will neglect this contribution to the total energy. 

D. Regularization of integrals and sums 

In calculating the vacuum contributions to the physical observables one has to deal with 
the appearance of divergent expressions. In this paper we consider two different regulariza- 
tion schemes, applying a momentum cutoff and using the Pauli-Villars regularization. 

In the first case we introduce a regulating function, which suppresses the contribution 
to integrals and sums at momenta k ^ A, where the scale A is determined by fitting the 
pion decay constant in free space (see Eq. ( p.l9|) . The regulating function we have chosen 
has the form 

Cik) = r„/V2WA21 1 - (2-38) 

exp [(A;2 - A2) /A^] + 1 ^ ^ 



In the Pauli-Villars regulating scheme p^ , on the other hand, the divergent contributions 
are eliminated through the subtraction 

K{x, x') -^ K{x, x') - K^'^ix, x'), (2.39) 

where K{x,x') is the propagator previously defined, while in i^^^(x,a;') the quark mass M 
has been substituted by the mass scale Mpv, obtained again by fitting the free space pion 
decay constant: 

/27rV2\ 



E. Pion decay constant 



Let's now consider the axial current. Its valence part is 



A^,{x)=^{x)^riMx). (2.41) 



We wish to calculate also the vacuum (i.e. from the Dirac sea at the one-quark-loop level) 
contribution to the axial current: This can be done by defining the generating functional 



W\a] 



dipdi/j exp i / d^x ip {ip - MU5) i/j - a'^^A^ 



(2.42) 



where a^ are classical axial sources coupled to the quantum axial current; then, the vacuum 
axial current can be obtained by means of a functional derivative with respect to the source 
as 



^a,vac\-^) 



5al{x] 



\aW\a\ 



<=o 



Calculation provides 

AN M^ r 

<vac(^) = ^^^ d'x'K{x,x') [d>^Mx')nx)-d'^rix')Mx)] . 

J TV 

By setting 

K{x, x') = Ko{r)6\x - x') + AK{x, x), 
where Ko{r) is given by Eq. ( p.27| ), one is able to write 

ANrM^ 



A^^{x)=^{x)—Y75Hx) + 



2 

ANcM^ 
P 

J IT 



J IT 



Ko{r) [d^^M^)rix) - d^rix)Mx)] 



d^x'AK{x,x') [9^00 (a;') 0" (a;) - 9^0''(x')0o(a;)] 



in which local and non-local contributions have been separated. 
Remembering now that the pion decay constant is defined as 

{0\A^,{x)\Ap)) = -W'fja,be-'^-'', 
one is able to obtain the in-medium pion decay constant as 

f'ir)=AN,M'Koir), 



(2.43) 



(2.44) 



(2.45) 



(2.46) 



(2.47) 



f2.48) 



that is expression (|2.31| ), which depends, in a medium, on the radial coordinate. 

An estimate of the average value of /^r at fixed density can be obtained by calculating the 
constant value of /tt.ws that would yield the same pion kinetic energy as the r-dependent 
one (compare Eqs. ' (|O0| ) and (glSOD ): 



iur 



lo'' dr f' (r) \r'9'\r) + sin' 9ir) 



J^dr r2^'2(r) + sin2 0(r 



(2.49) 
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F. Axial coupling constant 



The axial coupling constant, for a system with a finite pion mass is given by 

2^a = (pT \jdrAl\p]). 



(2.50) 



The matrix element between proton states requires one to adopt a collective quantization 
procedure in order to project out the correct quantum numbers. In the case of qa, — 
as discussed, e. g., in Ref. |Q, — it amounts simply to multiply the expression for the 
hedgehog axial current by the matrix element of the cranking operator, —1/3. Furthermore, 
as discussed in Ref. 1^^, the m.,^ = limit requires one to perform first the angular integral 
and then the radial integral in ( p.50| ), multiplying the result by a factor 3/2. 
Use of Eq. (|2^) yields 



9a 



R 



drr 



N, 



c / 2 



A-R 



J n 



Us'.^^ 



drAA 



3' 



(2.51) 



with 



drAA' 



IQnNrM^ 



Yl {/o("o)AfCo(ao, a'o)M<^o) + fi{ai)AKi{ai, a[)fi{a[)} 



a,a' 



and 



R 



/o(«o) = M'/^ I dr r' p^^{r) (cose(r) - 1) 

JO 

/i(ai) = M^/2 / dr r^ p^^{r) sin^(r). 

JO 

These are the expressions we shall employ in the next Section to evaluate the in-medium 
modification of the axial coupling constant. 



G. Equations of motion 

The equations of motion are found by minimizing the total energy ( |2.9| ) with respect to 
the Dirac fields, u and f , and to the chiral angle, 9. It is convenient to write them in a 
dimensionless form by letting q — > Mq and introducing 



X = Mr, u{x) = M'^/'^uir), v{x) = M~^/'^v{r). 
Minimization of etot = Etot/M = e^^i + elic then yields 



sin 6* M — (evai + cos 9)v 



du 

dx 

-— = — ( sin6')'i; + (cvai — cos H}u 

dx \x 



(2.52) 

(2.53a) 
(2.53b) 
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d'^e 2de sm29 N,.M'^ 



dx^ X dx 



+ 



x^ 

roo 



4vr/2 



-\ — sin 6 dqq 

71 Jo 

2 f^ 1 
cos 6 dqq 

-K JO 



K{q) 
K{0) 

K{q) ^ 
K{0) 



[(u^ — f ^) sin 9 + 2uv cos 6] 
3o{qx)C{q) 
ji{qx)S{q) = 0, 



where 



Ciq) 
S{q) 



dx{qx)^jo{qx)[cos9{x) — 1] 
dx{qx)'^ji{qx) sin6'(x), 



(2.53c) 

(2.54a) 
(2.54b) 



and K{q) is still given by (|2.22| ), but now k, k' are dimensionless and E^ ^* E^ = \/P~+T. 
This is a set of integro-differential equations that has to be solved iteratively as shown, for 
instance, in Ref. pT]] . 

The same procedure, applied to the system enclosed into the WS cell, yields 



-— = — sm 6' M — (evai + cos U)v 
dx 

— = -[ sm.6\v+ (cvai - cos 6)u 



(2.55a) 
(2.55b) 



dx"^ 



- + —lnKo{x) 



+ sm^— -- + cose 

Kq[x) Kq[x) 



dO Sm29 1 ., 2 ~2\ ■ n r.~~ m 

Tx-^^-^lMKjx)^^'' -.)sm^ + 2«.cos^] 

mx 



0, 



(2.55c) 



where Kq{x) is still given by (|2.27| ), but now Ea^/R — > -^a^/^ = J{ai/Xy + 1 and p, 



^ai\ 



Pai{x) = M '^^'^Kaji{aix/X), having set X = MR. In ( |2.55c|) we have also set 

AKi{ai,a[)fi{a[) (2.56a) 



Waix) = E 



-Pai [x] + — 1 

X dx 



W'^i^) = ^ E { ^^^^o(«0' «o)/o(«[,) + 2 



X dx 



AK2(a2,«2)/2(«2)|- 

(2.56b) 



H. Boundary conditions 

In the free space problem the boundary conditions on the fields are determined straight- 
forwardly. From inspection of the Dirac equations one has v{0) = and {u{x),v{x)} — > 0, 
X — i> cxd; finiteness of the energy requires 6{x) — >■ when x — * oo, whereas by choosing 
9{0) = n one fixes to unity the topological charge associated to the pion field. This is 
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sometimes interpreted as the baryon number, but not in the present model [Q, where it 
is connected to the number of valence levels that are pushed out of the Dirac sea. Baryon 
number is fixed by the normalization condition /q°° dxx'^p{x) = 1, where p{x) = u'^{x) + v'^{x) 
represents the (dimensionless) baryon density. 

If we consider now nuclear matter as a collection of hedgehog field configurations centered 
at lattice points, we have to impose periodic boundary conditions on the fields (Bloch's 
theorem). One then obtains a band structure, that is a continuous set of states and an 
energy gap above the highest energy state. The physical meaning of a band structure of 
the quark levels is not clear in the present context, since it implies long-range correlations 
among quarks and one should be aware that this (and similar) model does not account for 
confinement. We shall come back to this point later. 

In order to simplify the approach, the Wigner-Seitz approximation assumes that the cell 
be symmetric: Then, the state at the bottom of the band is also spherically symmetric and 
the wave functions are flat. To describe the other states, — and, in particular, the top of the 



band, — many different assumptions have been made in the literature (see, e. g., Ref. |T^ 
for a brief summary). In the calculations presented in Sect. |IT|, — besides the normalization 
condition /g dxx'^p{x) = 1, which fixes the baryon number, — we shall use three distinct 
sets of boundary conditions (sets I, II and III). 

In order to make contact with previous calculations, we follow for set I the choice of 



Ref. [19|, where the authors insist in maintaining unit topological charge inside the cell: 



^(0) = TT, 9{X) = 

v{0) = 0, v{X) = ("bottom" of the band) (2.57) 

u{X) = ("top" of the band), 

where v{X) = which implies u'{X) = 0. 

For set II we choose to impose "flatness" also on the chiral angle, as in Refs. PJT^. Since 



in general one has 0{X) 7^ at the boundary, from inspection of the Dirac equations one 
sees that v{X) = no longer implies u'{X) = 0. We have chosen to impose the physically 
motivated constraint 

p\X) = 0, (2.58) 

that is we require the flatness of the baryon density at the cell boundary. Set II then turns 
out to be 

^(0) = TT, e'{X) = 

Xcos^(X) T Jx^-2Xsme(X) 

m = 0, v{X) = ' .V '—u{X) ("bottom"/"top"). (2.59) 

X smt'(A) — 2 

For set III we require the flatness of the baryon density together with the requirement of 
unit topological charge: 

^(0) = TT, e{x) = 

m) = /^W = ("bottom") 

v[^) u, |^(x) = -Xu(X) ("top") ■ ^^■^^> 
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One might be tempted to interpret the two solutions of the equation ( p.58| ) as corresponding 
to the bottom and top of the energy band. Indeed, from set III one sees that v{X) = 
corresponds to the bottom of the band in set I, whereas in the non-relativistic limit {v — *• 0) 
the second condition in ( p. 601 ), that is —Xu{X) = v{X), reduces to u{X) = as in set I. 



A word of caution is necessary in analyzing the boundary conditions. In our view one 
cannot accept without question the presence of a band of quark states. The presence of a 
band, in fact, would be affected by confinement, which is absent in the present model. In 
free space, where the quarks are deeply bound in the ground state of the chiral fields, this 
shortcoming is not crucial (of course, the study of the highly excited states of the nucleon 
would then be problematic). In the medium, however, because of the lack of confinement, 
quarks exhibit unrealistic long-range correlations. As a result, one observes a relatively 
large probability of having a quark sitting at the surface of the WS cell. 

In the next Section we shall see that when the density of the medium increases the quark 
density tends to be more concentrated in the interior of the bag for the lowest end of the 
energy band (lowest when all effects have been included), whereas the opposite happens on 
the upper end. As a consequence, quarks sitting in the upper part of the energy band would 
be more affected by the confining forces than the ones in the lower part. Since the confining 
forces would tend to reduce the quark density at the boundary, the net result would then 
be mainly a lowering, — and hence a narrowing, — of the highest end of the band. 

Another problem one should cope with in identifying the top and bottom energy levels is 
posed by the presence of the spurious center of mass energy. Although this contribution has 
been neglected in most existing calculations of in-medium properties in the Wigner-Seitz 
approximation, it turns out not only to be sizeable but also to affect the relative position of 
the "top" and "bottom" levels, as it will be discussed in the next section. 

In the following, for definiteness, we shall follow the nomenclature adopted in the lit- 
erature and we shall label "bottom" and "top" the two solutions of each set of boundary 
conditions, as indicated in ( ^.57| ), ( ^.58| ) and ( |2.60| ). 

III. RESULTS 



The numerical results have been obtained by integrating the equations of motion (|2.55|) 



for the quarks and the meson fields, using an iterative procedure as, e.g., in Ref. 1^^ and 
moving from larger to smaller values of R. For a given value of R, i.e. for a given density, 
a self-consistent solution has been found by using as initial ansatz the self-consistent chiral 
profile obtained at the previous value of R. The non-local term in Eq. (|2.55c| ) has been 
switched on adiabatically in order to allow a better convergence^. 

At the smallest density, corresponding to i? = 5 fm, an exponential profile 6{r) = 
n exp{~r/ro) has been used as initial ansatz. A dynamical (constituent) quark mass 
M = 350 MeV has been assumed in the calculations. This is a value suggested by the phe- 



^The sums over the modes in the orthonormal and complete basis inside the spherical cell have 
been restricted to / < 15 and to the first 30 roots a/, which provides a good degree of accuracy. 
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TABLE I. Comparison between the energies (in MeV) obtained in free space and in the WS 
cell using the boundary conditions of sets I-III for R = 5 fm. 



Eval 



Ekh 






(2) 



i?kin + E(') 



E, 



tot 



free 

II 

I and III 



452.4 
443.6 
452.9 



922.4 
921.8 
917.5 



-387.6 
-386.2 
-389.6 



534.8 
535.6 
527.9 



987.2 
979.2 
980.8 



nomenology of the single nucleon 0-0], which also turns out to be in the range (300 — 400 
MeV) where the second order expansion for the effective action works well p7[| . 

The divergences that appear both in free space (in the momentum integrals) and in 
the medium (in the sums) have been regulated as explained in Sect. [II D| , using both a 
regulating function and the Pauli-Villars regularizations, whose parameters have been fixed 
by fitting the free space value of the pion decay constant. For M = 350 MeV, the cutoff 
in the regulating function turns out to be A = 500 MeV, whereas the mass scale Mpv of 
the Pauli-Villars regularization is given by Eq. ( 2.40| ). The two regularization schemes yield 
qualitatively similar results and in the following we shall display only the outcome from the 
regulating function approach. 

Before discussing the results a comment on the fiat basis we have adopted is in order. 
In fact, this basis contains a zero momentum state, that is, a term constant in space and 
proportional to R^^^"^. When i? — > 0, it would give rise to divergences. As explained 
in Appendix ^, the appearance of such a mode is peculiar of the flat basis: One might 
introduce bases infinitesimally close to the flat one, in which the zero mode is absent and its 
strength is distributed among all the other modes. In these alternative bases the effect of 
the zero mode would be given by an infinite sum of infinitesimal contributions; since we are 
regularizing the sums, only a finite number of modes enters into the calculation of physical 
quantities and the contribution to them of the redistributed strength of the zero mode is 
infinitesimal. Hence, we can cure the divergence simply by dropping the zero mode. 

In order to exhibit the convergence of the WS calculations to the free space results, let us 
start by comparing, in Fig. |I|, the chiral angle obtained by solving the equations of motion 
in free space (solid lines) to the WS solutions corresponding to i? = 5 fm and obeying the 
boundary conditions of set II (dashed line) and of sets I and III (dotted line). Note that at 
such a low density, the top and bottom energy levels practically coincide. In both cases a 
regulating function has been used to regulate sums and integrals. As expected, the difference 
between the solutions inside the WS cell and the one in free space is barely noticeable at 
this density. By looking at Table |I| we also notice that the free space energies are recovered 
with good precision, using any of the boundary conditions. 

An important difference between our work and previous calculations employing chiral 



quark models p!7|-pO| is due to the inclusion in our calculations of non-local effects stem- 
ming from the vacuum contribution (at the one-quark- loop level). In Fig. § we display, as 
a function of R, the WS cell total energy (solid line), separated in the valence (dotted line) 
and vacuum (dashed c line) terms, for the "bottom" solution of set I. Also shown are the 
separated local, — that is, kinetic, — (dashed a line) and non-local (dashed b line) contribu- 
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FIG. 1. The chiral angle in free space (solid line) and in the WS approximation at i? = 5 fm, 
using the boundary conditions of set II (dashed) and of sets I and III (dotted). 
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FIG. 2. Total energy of the WS bag (solid) for the "bottom" solution of set I. Also shown are 
the valence (dotted) and vacuum (dashed c) contributions; of the latter one, we display also the 
local (dashed a) and non-local (dashed b) components. 
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FIG. 3. Total energy of the WS bag. The left panel corresponds to the boundary conditions 
of set II (solid and dashed lines for the "top" and "bottom" solutions, respectively), whereas the 
right panel corresponds to the boundary conditions for the "bottom" solutions of sets I and III 
(dashed) and for the "top" solutions of sets I (dotted) and III (solid). 
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FIG. 4. Binding energy obtained by taking out the spurious center-of-mass energy contribution 
and subtracting the energy at the lowest density (here R = 5 fm). The left panel corresponds to 
the boundary conditions of set II (solid and dashed lines for the "top" and "bottom" solutions, 
respectively), whereas the right panel corresponds to the boundary conditions for the "bottom" 
solutions of sets I and III (dashed) and for the "top" solutions of sets I (dotted) and III (solid). 
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tions. The local term displays a behavior similar to the results of, e. g., Ref. pO[: Actually, 



in the chiral quark model employed in that paper, the kinetic meson contribution is present 
at the classical level, whereas in our case it is dynamically generated from the vacuum. 
This implies, as we saw in Sect. ||, a r-dependent pion decay constant (see Eqs. (|2.30|) and 
( p.31| )); We have however checked that setting /tt^ws constant in our calculation, the results 
of Ref. [^ are recovered. On the other hand, the non-local vacuum contribution provides 
substantial attraction and displays a moderate dependence on R; however, it turns out that 
the valence and kinetic meson terms compensate each other to a large extent yielding a total 
contribution without any minimum and the non-local term is then instrumental in order to 
get the (rather shallow) minimum displayed by the solid line. 

By looking at Fig. ^ one notices that solutions to the equations of motion are no longer 
found below R ^ 1.4 fm. This depends also on the choice of the boundary conditions and 
in the other cases discussed below we shall see that solutions are found till R ^ 1.1 fm. 
In Ref. |2^ solutions have been found till much higher densities {R ~ 0.4 fm). In the 
case of Fig. 2 in that paper, — where a simple Lagrangian containing only terms up to 
second order in the pion field is used, — this is due to the larger value for the constituent 
quark mass M chosen in that work (in their notation M = gf.,^ ^ 550 MeV, where g is 
the quark-meson coupling constant). We have checked that by dropping the non-local term 
and by increasing M their results can be recovered. Indeed, from Eq. ( |2.55| ) one sees that 
the equations of motion depend on the combination X = MR: By increasing M , one can 
lower the minimum value of R. However, the authors of Ref. pO| are able to find solutions 



at higher densities by using more complex Lagrangians containing terms of higher order in 
the pion field. Following that path in our model would imply going beyond the two-point 



approximation to the effective action of Sect. II B 



In Fig. 1^ we display the WS cell total energy as a function of the cell radius, i.e. of the 
density. The boundary conditions of set II and of sets I and III have been used in the left 
and right panels, respectively. In Fig. ^ we notice the presence of a very shallow minimum 
around a density corresponding to i? ~ 2 fm; this minimum is deeper for the "bottom" 
solutions. 

The occurrence of saturation in nuclear matter cannot however be stated by simply 
looking at these figures, because here the spurious energy contribution due to the center- 
of-mass motion has been neglected. This, as a matter of fact, is a well known problem 
associated to the mean field approximation. An estimate of this effect in the chiral quark 
soliton model has been obtained in [^, including also the vacuum (mesonic) contributions 



to the center of mass motion. The findings in that paper, — that valence terms dominate 
as long as E^ai > 0, — make us feel confident in retaining only the latter in our estimate. 



It reads 

EcM 



{P') 



2Eu 



N f /-^ 



2K 



tot 



R^p'{R) - J dr [r\u'\r) + v'\rf) + 2v\r)] I . (3.1 



Of course, the validity of this assumption at finite density can only be checked through an 
explicit calculation of the vacuum terms, which is however beyond the scope of the present 
analysis. 
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FIG. 5. Dependence of the baryon density upon the cell radius for the "top" (solid) and 
"bottom" (dashed) solutions of set II. 
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FIG. 6. Decomposition of the total energy (solid) into valence (dot) and vacuum (dashed) 
contributions, for the "top" solutions of set II (left panel) and III (right panel). 
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In Fig. ^ we plot the binding energy for the system taking out the spurious contributions 
stemming from the motion of the center of mass. In order to minimize the numerical 
uncertainty, the binding energy has been obtained by subtracting the total energy at the 
lowest considered density {R = 5 fm). Interestingly, we find that, when the center-of-mass 
motion is taken out, a stronger minimum in the energy is found, roughly at the same density 
as in Fig. ^, namely i? ^ 1.8 fm. Moreover, the boundary conditions of sets I and III provide 
more binding that those of set II. Note that in nuclear matter one should have a binding 
energy of about —16 MeV at i? = 1.1 fm. 

We also notice that the "top" solutions of the various sets of boundary conditions provide 
more binding than the corresponding "bottom" solutions. The reason is easily understood 
by looking, for example, at Fig. |[ 

Here we plot the value of the baryon density at the surface of the cell, i.e. v?{R) +f ^(i?), 
as a function of the cell radius itself, for the boundary conditions of set II (the solid and 
dashed lines corresponding to the top and bottom solutions respectively); sets I and III 
show a similar behavior. Since the top solution corresponds to a configuration in which 
the quarks are more "compressed" inside the cell, a larger kinetic energy, — and therefore 
a larger center-of-mass motion, — is associated with it. This problem is of course absent 
in solid state physics, the original field of application of the Wigner-Seitz approximation, 
because the electron mass is indeed completely negligible with respect to the mass of the 
ions, which form the periodic structure. In the present case, even admitting the existence 
of a periodic structure in nuclear matter, the center-of-mass motion would be neglegible 
only in the large Nc limit {N^ -^ oo), given the dependence of the total energy and of the 
center-of-mass energy on the number of colors, as 0{Nc) and 0(iV°) respectively. On the 
other hand, for Nc = 3, the center-of-mass energy will vary, for each solution inside the 
band, by an amount comparable with the width of the band itself. Hence, the calculation 
of a reliable band structure within this model is more delicate. 

In Fig. ^ the total energy corresponding to the "top" solutions of set II (left) and HI 
(right) is plotted, and decomposed into the valence (dotted) and vacuum (dashed) contri- 
butions. We observe that, at very low densities, the two components bear approximately 
the same strength, whereas at larger densities, below R = 2 fm, the valence contribution 
becomes dominant. It is also interesting to display, as a function of R, the ratio between the 
valence contributions that come from the quark kinetic term in the Lagrangian and from the 
quark coupling to the mean field. This is done in Fig. 0: As expected, this ratio increases 
with the density. 

In Fig. Ijthe average value of the pion decay constant in the cell, defined in Eq. ( |2.49|) , is 
plotted as a function of the cell radius, for the "top" solution of set II; for all the other sets of 
solutions the behavior is very similar. The pion decay constant decreases by increasing the 
density, going in the direction of a partial restoration of chiral symmetry, although the lack 
of solutions beyond roughly the standard saturation density of nuclear matter, prevents one, 
— at the present stage of development of the model, — from drawing firmer conclusions. 
A reduction of /^r in matter is found both in linear sigma and Nambu-Jona-Lasinio models 
(see, e. g., Ref. |^ for a list of references). 



In Fig. ^ the isoscalar mean radius is plotted as a function of the cell radius, for the 
"top" (solid) and "bottom" (dashed) solutions of set II; the other sets of solutions give 
very similar results. We observe that this quantity is extremely sensitive to the choice 
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FIG. 7. The ratio between the kinetic and potential components of the valence energy as a 
function of the cell radius, using the "top" boundary conditions of set II (solid) and set III (dashed). 
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FIG. 8. Dependence of the pion decay constant on the the cell radius for the "top" solution 
of set II. 
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FIG. 9. Isoscalar mean square radius, as a function of the cell radius. The solid line cor- 
responds to the "top" solution of set II, whereas the dashed line corresponds to the "bottom" 
solution of Set III. 
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FIG. 10. Chiral angle as a function of the distance from the center of the soliton for R = 5 
(solid), 2 (dashed) and 1.1 fm (dotted), respectively. The left panel corresponds to the "top" 
solution of set II, whereas the right one corresponds to the "top" solution of set III. 
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FIG. 11. Large and small components of the Dirac spinor as a function of the distance from 
the center of the soliton for R = 5 (solid), 2 (dashed) and 1.1 fm (dotted), respectively. The left 
panel corresponds to the "top" solution of set II, whereas the right one corresponds to the "top" 
solution of set III. 
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FIG. 12. Baryon density as a function of the distance from the center of the soliton for R = 5 
(solid), 2 (dashed) and 1.1 fm (dotted), respectively. The left panel corresponds to the "top" 
solution of set II, whereas the right one corresponds to the "top" solution of set III. 
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FIG. 13. The axial coupling constant as a function of R, for the "top" boundary conditions of 
set II (left panel) and set III (right panel). The valence (dashed), local mesonic (dotted), non-local 
mesonic (dot-dashed) and total (solid) contributions are displayed. 



of "top" or "bottom" boundary conditions, the reason being the same as aheady noticed 
when discussing Fig. ^ The "bottom" solution in fact corresponds to a configuration in 
which the quarks are more loosely packed inside the cell. One could in principle define a 
mean radius averaged over the solutions distributed in the band, assuming, for instance, 
a uniform distribution. From Fig. |^ one would expect such an averaged radius to have a 
moderate density dependence. However, as already mentioned, we expect confinement to 
affect the band width, and expecially the "bottom" border (actually, the upper end), thus 
unbalancing the distribution in favor of the lower curve of Fig. ^ and yielding a shrinking of 
the soliton. This would be at variance with findings of other models. 

For completeness we display in Figs. [TO, 111] and O, as a function of the distance from 



the center of the WS cell and for three values of R, the chiral angle, the large and the 
small components of the Dirac spinor and the baryon density, respectively, using the "top" 
solutions of sets II and III. Similar results apply for the other solutions. 

A quantity which is more sensitive to the choice of boundary conditions is qa, as one can 
see in Fig. |13|. The local mesonic component varies noticeably in the different cases, even 
at large R. In Table || we display the value of the different components oi qa 9.t R = 5 fm, 
compared to the results for the free soliton caseQ. One clearly sees that the quark and 



^ As noted by the authors of Ref. |27|, the two-point approximation to the effective action works 
for QA to a lesser extent than for the energy and the apparent agreement with the experimental 
value should be regarded as accidental, since the self-consistent calculation (see, e. g., Ref. 0) 
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TABLE II. 








quark 


local mesonic 


non-local mesonic 


total 


free 


1.07 


0.50 


-0.35 


1.22 


II 


1.06 


0.76 


-0.34 


1.48 


III 


1.06 


1.03 


-0.34 


1.76 



TABLE III. The axial coupling constant qa in free space and in the Wigner-Seitz approxima- 
tion at i? = 5 fm, for the "top" boundary conditions of sets II and III. 

non-local mesonic components have already reached the asymptotic value, whereas the local 
mesonic term is still higher. The reason for this behavior can be traced back to the slow 
decay of the chiral angle in the chiral limit {0{r) ^ 1/t^), which allows sizeable contributions 
from large distances. Although in principle one could continue the WS calculation at larger 
radii, one would then need to increase the number of states included in the orthonormal basis 
employed in the expansion. However, it is likely that, when chiral symmetry is explicitly 
broken, the exponential decay of the chiral angle would grant a faster convergence and 
suppress the sensitivity to the boundary conditions. Increasing the density, the mesonic 
contributions to qa, both local and non-local, rapidly decrease and most of the strength is 
now carried by the quarks alone. 

IV. OUTLOOK AND PERSPECTIVES 

In this paper we have applied the Wigner-Seitz approximation to the chiral quark-soliton 
model of the nucleon. This model, complemented by the WS approximation, provides a sim- 
ple yet interesting framework for studying possible modifications of the nucleon properties in 
nuclear matter. It is important to remark that we are actually dealing with a parameter-free 
model, since the only free parameters (the constituent quark mass M and the regularization 
scale) have been fixed in free space. Note that M could, in principle, be calculated dynam- 
ically, as done, e. g., in the Nambu-Jona-Lasinio model [^; this would require the solution 
of a gap equation and would probably lead to an in-medium suppression of M, similar to 
what is found for fj^. 

The same method of accounting for in-medium effects has already been applied to a 



number of microscopic models of the nucleon ||T3|-p3|. With respect to previous analyses, we 



have for the first time consistently calculated effects stemming from the Dirac sea, i. e., from 
excitations of virtual quark-antiquark pairs, including their dependence upon the density of 
the system. 

Vacuum fluctuations manifest themselves in two ways. Firstly, by giving rise to a pion 



gives a smaller value. 

^By increasing the density, only virtual states of higher energy can be excited and, eventually, in 
a very dense system, only states lying well above the typical momentum cutoff would be available: 
The quark condensate, (qq), is thus expected to decrease at finite densities. 
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decay constant dependent on the distance from the center of the bag, which in turn modifies 
the pion kinetic contribution with respect to chiral models where the pion is exphcitly 
included at the classical level p!8|- p0| , p!7| (see Eqs. ( p.30| ) and ( |2.31| )): The contribution to 
the energy from this term is however qualitatively similar to the one of previous analyses 
p!8|- p0| , p!7| . Secondly, vacuum fluctuations generate a new non-local, attractive contribution, 
whose density dependence, albeit moderate, is however relevant in order to bind the system, 
given the compensation one observes with increasing density between the valence quark and 
kinetic pion contributions. 

Another effect that has usually been neglected in previous calculations is the spurious 
center-of-mass motion: We have found it to be important not only to give, of course, a more 
realistic estimate of the energy, but also to determine the relative position of the top and 
bottom ends of the quark energy band in the medium. These two levels correspond to specific 
boundary conditions and we have found the role of the boundary conditions generating the 
upper and lower levels to be exchanged with respect to previous calculations. Note that the 
intersection of an occupied band with an empty one is often interpreted as the onset of color 
superconductivity. 

We have also explored the effects stemming from different choices of boundary conditions. 
All the cases discussed above display a similar qualitative behavior, although there are 
definitely quantitative differences, especially for the axial coupling constant, which is however 
very much affected by the long-range behavior of the chiral field in the (massless pion) chiral 
limit. Also the two regularization schemes we have adopted turned out to give qualitatively 
similar results. 

We have calculated a few physical quantities, such as f-,, and (r^)/=o- The pion decay 
constant has been found to decrease with increasing density, pointing to a partial restoration 
of chiral symmetry in the medium. The isoscalar mean square radius, on the other hand, has 
been found to depend heavily on the position in the energy band, leaving open the question 
whether an average nucleon will swell or shrink, given the present uncertainties in assessing 
the band structure. 



Unlike previous calculations with chiral models [p]8 H20|jr^ , we have been able to find bind- 
ing, but at a density much lower than the standard saturation density of nuclear matter; 
moreover, solutions of the equations of motion disappear roughly below the latter density. 



The authors of Ref. I^O] have been able to find solutions at high densities, — still keeping 
realistic values for the constituent quark mass parameter M , — by incorporating in their 
chiral Lagrangian terms of higher order in the pion field. In the present model this would 
correspond to dropping the two-point approximation to the effective action in evaluating 
vacuum fluctuations. This is probably the most needed development of the present calcula- 
tion, — employing, e. g., the numerical algorithm of Ref. [^, — since only in the free case 
it has been explicitly verified that the two-point approximation works reasonably well [^ 



Finally, another feature that is lacking at the present stage is the projection of good 
spin-isospin quantum numbers out of the hedgehog state. Since the momentum of inertia is, 
in general, density-dependent, the projection is likely to affect the position of the minimum 
in the energy. 
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APPENDIX A: THE PROPAGATOR K{Q) 



We derive here the expression for the propagator K{q), introduced in Sect. p.lB| . As a 
first step we perform the angular integrals in Eq. ( p.l8| ), getting 



K{q) 



1 foo roo 



271^ Jo 



dkk' 



,,,, n C[ik + k')/2] 



dk'k 



EkEk' Ek + Ek' 



drr^jo{qr)jo{kr)jo{k'r). (Al) 



The inner integral can be done analytically ||3^, yielding 

3/2 I 



J drr^jo{qr)joikr)jo{k'r) = i-j r-rri J dr ^/rJl/2{qr)Jl/2{kr)Jl/2{k'7 
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where the function A(a;, y, z) vanishes whenever it is not possible to build a triangle of sides 
x, y and z, and is equal to 1 otherwise. 
One then obtains 



K{q) 



iTi'^q Jo 



Ek J\k-q\ Ek' Ek + Ey 



(A3) 



APPENDIX B: A BASIS FOR FLAT FUNCTIONS 

In free space {R — >■ cxd) an orthonormal and complete set of states can be chosen as 

^fcoMm(x) = -e'^'^^Vrmjiikr). (Bl) 

TT 

We want to build an orthonormal and complete basis inside a sphere of radius i?, having 
the same form of ( pi|) . While in that case the momentum could take a continuum set of 
values, now only a discrete set of momenta will be allowed. 

We start considering the equations for the spherical Bessel functions, corresponding to 
two different momenta, fci = a/R and /c2 = /3/-R: 

= 
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(B2) 
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Multiplying both expressions by ji{ar/R) and ji{f3r/R), respectively, integrating over r 
between and R, taking the difference of the two equations and, finally, integrating by 
parts, one gets 

- I^"" {a' - (3') ^j, (^) J, (^^) = R [ajl{a)MP) - mf3)Ma)] (B3) 

The above equation states the orthonormality of the elements of the basis, provided that 
the following condition is met: 

ajl{a)ji{P)-f3jl{P)jiia) = 0, a^(3. (B4) 

The choice of the boundary conditions fulfilled by the elements of the basis will therefore be 
constrained by this requirement. It is easy to convince oneself that three different possibilities 
are available: 

3i{a) = (B5a) 

fiia) = (B5b) 

"Ji(") = VJi{a), (B5c) 

where r] will is a constant parameter. In the limit rj = and rj ^ oo the first two cases are 
recovered, respectively. 

Notice that, in the limit R -^ oo, — where the modes form a continuum, — all these 
boundary conditions become equivalent. 

We then introduce the normalized functions 

f \ - ■ f^i^ 



2 _A 



c^ifiic^i? - Jiic^i) (-2«/j;(az) + (/(/ + 1) - afMai))] ' , (B6) 



such that 

drr^paiir)pf3Xr)=6a,f3r (B7) 

The completeness of the basis allows one to write the following representation for the Dirac 
delta function inside the WS cell: 

EP".WP.,(r') = ^^^7^- (B8) 

It is simple to derive an expression for the Green's function of the operator (□ + M^) in the 
WS cell. It reads 



Gix,x') = Y.jy-^Y: Ynmrme^''^'^"'''^ avl^'+M^'- p ■ ^^^^ 



The fiat basis is obtained, as we said above, in the limit r^ -^ 0; in this case, however, the 
lowest energy mode, — corresponding to zero momentum for / = 0, can only be obtained for 
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free space Wigner-Seitz boundary conditions 

k — > OLi/R 

jijkr) —2_ Pai {r) 

TABLE IV. Rules for going from free space to the Wigner-Seitz cell. 



a strictly vanishing 77 (in fact, it corresponds to a negative root of Eq. ( [B5c| ) at finite values 
of 7]). This means that the strength of this mode has to be redistributed, for finite 77, over 
all the other modes. Let us call (3 and (3 the modes such that 

j;(/5)=0, M~P)=r^jm- (BIO) 

There is a one-to-one correspondence between the modes in the two bases, in such a way 
that for each /3 7^ one can define (3 = (3 — e^, where e/3 <^ 1 if 77 ^ 1. 
Writing down the transformation from one basis to the other, namely 

P/3W = IIc^^P/3W. (Bll) 

with some algebra one can verify that, for / = and /? = 0, one has 

J 



co^ = -^^P, (B12) 



which indeed proves that the zero mode is now redistributed over all the other modes, in 
the quasi-flat basis. 

Calculations in the two bases are equivalent (to order 77) only when considering convergent 
quantities. Quantities that need regularization are therefore different in the two bases, since 
only a finite number of modes will contribute. The fiat basis is the one we have chosen in the 
calculations of this paper, since it is the most convenient in order to impose the physically 
motivated boundary conditions discussed in the text. The zero mode, however, gives rise to 
divergences when R ^ and one could in principle eliminate this problem by employing a 
quasi-flat basis infinitesimally close to the flat one. On the other hand, the very fact that 
for regularized quantities only a finite number of modes is relevant, allows us to keep the 
flat basis and simply discard the zero mode, since the error introduced in this way will be 
infinitesimally small. 



Finally, we display in Table IV a few simple rules that allow one to rewrite in the WS 



basis quantities expressed in the free basis (p31| 
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